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Abstract 



We study pair-production as well as the triple self-couplings of the neutral Higgs bosons of 
' the Minimal Supersymmetric Standard Model (MSSM) at the Future International Linear e'^e" 

Collider (ILC) and Compact Linear Collider (CLIC). The analysis is based on the reactions e~^e~ — 
on : m.h., tth.h, with = h, H, A. We evaluate the total cross-section for both bbh.h., tih^h. and 

' calculate the total number of events considering the complete set of Feynman diagrams at tree- 

' level. We vary the triple couplings nXhhh, i^^Hhh-, i^^hAA-, i^^haa-, i^^hhh and hXhhh within the 

o : 

■ range k = —1 and +2. The numerical computation is done for the energies expected at the ILC 

with a center-of-mass energy 500, 1000, 1600 GeV and a luminosity 1000 fb-^. The channels 

^ ' e+e~ — > bbhihi and e~^e~ tthihi are also discussed to a center-of-mass energy of 3 TeV and 

H \ 

luminosities of 1000 fb~^ and 5000 fb^^. 

PACS numbers: 13.85.Lg, 12.60.Jv 
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I. INTRODUCTION 



The search for Higgs bosons is one of the principal objectives o 



energy coUiders, such as the International Linear Collider (ILC) [l 



Ref. 



Collider (CLIC) |5( and Large Hadron Collider (LHC) 



■present and future high- 



m 



Compact Linear 



It has been demonstrated in 



that physics at the LHC and e^e~ ILC will be complementary to each other in 
many respects. In man y c ases, the ILC can significantly improve the LHC measurements. 
The Higgs boson 



lanycc 
0,0, 
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13| plays an important role in the Standard Model (SM) 



Iq because it is responsible for generating the masses of all the elementary particles 



(leptons, quarks, and gauge bosons). However, the Higgs-boson sector is the least tested in 
the SM, in particular the Higgs boson self-interaction. If Higgs bosons are responsible for 
breaking the symmetry from SU{2)l x U{1)y to U{1)em, it is natural to expect that other 
Higgs bosons are also involved in breaking other symmetries. One of the more attractive 



extensions of the SM is Supersymmetry (SUSY) |l7l.ll8l. Il9|. mainly because of its capacity to 
solve the naturalness and hierarchy problems while maintaining the Higgs bosons elementary. 

The theoretical framework of this paper is the Minimal Supersymmetric extension of the 
Standard Model (MSSM), which doubles the spectrum of particles of the SM, and the new 



free parameters obey simple relations. The scalar sector of the MSSM [20|, |2l| requires two 
Higgs doublets, thus the remaining scalar spectrum contains the following physical states: 
two CP-even Higgs scalares (/i°,iy°) with < Mh, one CP-odd Higgs scalar (y4°) and a 
charged Higgs pair (H^). The Higgs sector is specified at tree- level by fixing two parameters 
which can be chosen as the mass of the pseudoscalar and the ratio of vacuum expectation 
values of the two doublets tan/5 = V2/V1; then, the masses M^, Mh and Mh± and the 
mixing angle of the neutral Higgs sector a can be fixed. However, since radiative corrections 
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29|, it 



produce substantial effects on the predictions of the model [22, l23|, 
is necessary to also specify the squark masses, which are assumed to be degenerated. 

In particular, all the triple self-couplings of the physical Higgs particles can be predicted 
(at the tree level) in terms of Ma and tan/3. Once a light Higgs boson is discovered, the 
measurement of these triple couplings can be used to reconstruct the Higgs potential of the 
MSSM. 

The triple Higgs self-couplings can be measured directly in pair-production of Higgs 
particles at hadron and high energy e^e~ colliders. In proton collisions at the LHC [7|, 
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Higgs pairs can be produced through double Higgs-strahlung off W and Z bosons, WW and 
ZZ fusion, and gluon-gluon fusion. The triple Higgs boson couplin gs a re involved in a large 
number of processes at e+e" collider [sol, H, Q, H, Q, Q, Q, HH, Q, 0,0,13: 



double Higgs-strahlung : e~^e~ — > ZHiHj and [ifjj = h,H], 

triple Higgs production : e~^e~ — ^ AHiHj and AAA, (1) 
lyVT fusion : e~^e~ — ^ Ve^eHiHj and UgUgAA, 

these three-body processes have been evaluated extensively. However, the inclusion of four- 
body processes with heavy fermions b and t {e~^e~ —>■ bbhihi, tihihi, hi = h, H, A) is important 
in order to know its impact on three-body channels and also to search for new relations that 



the 
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may have a clear signature of the Higgs boson production. In the studied proces ses, 
MSSM Higgs bosons are radiated by b(b) and t(t) quarks at future e"'"e~ colliders [43, 

481] with a cm. energy in the range of 500 to 1600 GeV, as in the case of the 
ILC 111 and of 3 TeV to the CLIC machine |5|. 

In this paper, we study the pair-production as well as the triple self-couplings of the Higgs 




le Future International 



bosons of the Minimal Supersymmetric Standard Model (MSSM) at t 

Linear e~^e~ Collider (ILC) 1] and Compact Linear Collider (CLIC) jsj. The analysis is based 
on the reactions e~^e^ bbhihi, tihihi with hi = h, H, A. We evaluate the total cross-section 
for both bbhihi, tihihi and calculate the total number of events considering the complete 
set of Feynman diagrams at tree- level. We vary the triple couplings i^Xhhh, f^^uhh, i^^hAA, 
i^^HAA, i^^hHH and hXhhh within the range k = — 1 and +2. The numerical computation is 
done for the energies expected to be available at the ILC with a center-of-mass energy 500, 
1000, 1600 GeV and a luminosity 1000 fb~^. Our analysis is also extended to a center-of- 
mass energy 3 TeV and luminosities of 1000 fb~^ and 5000 fb~^. We consider tan/5 = 35 
and Ma = 400 GeV. 

The Higgs couplings with quarks, the largest couplings in the MSSM, are directly ac- 
cessible in the processes where the Higgs boson is radiated off bottom quarks, e~^e~ — >■ 
bbhh, bbHH, bbAA as well as in the processes e^e~ — tihh, tiHH, tiAA where the Higgs 
boson is radiated off top quarks. These processes depend on the Higgs boson triple self- 
couplings, which could lead us to obtain the first non-trivial information on the Higgs 
potential. We are interested in finding regions that could allow the observation of the 



bbhh, bbHH, bbAA and tthh, ttHH, ttAA processes at future linear e~^e colliders energies: 
ILC and CLIC. We consider the complete set of Feynman diagrams at tree-level (Figs. [UlS]) 



and use the CALCHEP 



49| packages to evaluate the amphtudes and cross-sections of the 



processes e~^e~ bbhihi,tthihi. 

This paper is organized as follows: In Sec. II, the Higgs boson interaction Lagrangian and 
the self-couplings are presented. In Sec. Ill, we study the triple Higgs boson self-coupling 
through the processes e+e~ —>■ bbhh{HH, AA) and e+e~ — > tthh{HH, AA) at future e+e~ 
colliders energies (ILC/CLIC) and, finally, we summarize our results in Sec. IV. 



II. HIGGS BOSON INTERACTION LAGRANGIAN 

The Higgs sector of MSSM includes five physical fields: two neutral CP-even Higgs scalar 
{h^,H^), one neutral CP-odd Higgs scalar (A^) and a charged Higgs pair (H^). The Higgs 
boson interaction lagrangian has the form 

r-Higgs _ ^(3) ^(4) /„n 
'-Int ~ '-Int "T '-Int^ \^ ) 

where df^^ is the lagrangian of the triple Higgs boson interactions, while Cf^^ is the la- 
grangian of the quartic Higgs boson interactions. Explicitly, CjJ^^ is: 



£(3) 



A 



hhh 



Int 



3! 
2! 



hhh + ^hhH + ^hHH + ^HHH + ^hAA 



2! 2! 3! 

HAA + XhH+H-hH+H^ + Xhh+h-HH+H~. 



2! 



(3) 



The triple Hig gs s elf-coup 
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ing of the MSSM can be par ameterized as 
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^hhh 

^Hhh 
XuHh 

Xhhh 



3 COS 2a sin(/5 -|- a) + 3 



e cos a 
M| sin (3 



cos a, 



e sm a 



2 sin 2a sin(/3 + a) — cos 2a cos(/5 + a) + 3-;^^^ — cos^ a, 



Mi sin p 



—2 sin 2a cos(/5 + a) — cos 2a sin(/5 + a) + 3 



e cos a 
M| sin (3 



sm a, 



3 cos 2a cos(/3 + a) + 3 



e sm a 



Ml sin p 



sm a, 



(4) 
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A/xAA = COS 2/3 sm(/3 + a) + . cos /:/, 

M| sm p 

. \ e sina 2 ^ 

AffAyi = -cos2^cos(/3 + a) + — cos P, 

Mz sm p 

where e ~ ^o^\ ^ log ^ and the connection between the mixing angles a and Q is ffiven 
by 



tan 2a = tan 2/5— 5 ^ with - - < a < 0, (5) 

'^M2 -M| + e/ cos /3 2" " ' 



function of Ma and tan /5. 



III. HIGGS BOSONS PAIR PRODUCTION WITH MODIFIED TRIPLE SELF- 
COUPLING 

In this section we present numerical results for e^e" bbhh{HH, AA) and e~^e~ — *• 
tihh{HH, AA) with double Higgs bosons production. We carry out the calculations using 
the framework of the Minimal Supersymmetric Standard Model at the future linear e'^e~ 



colliders. We use the CALCHEP 49|] packages for calculations of the matrix elements and 
cross-sections. These packages provide automatic computation of the cross-sections and 
distributions in the MSSM as well as their extensions at tree-level. We consider the high 
energy stage of a possible Future Linear e~^e~ Collider (ILC) with ^/s = 500, 1000, 1600 
GeV and design luminosity 1000 fb^^. For the numerical computation, we have adopted 
the following parameters: the angle of Weinber sin^ 9\y = 0.232, the mass (m^ = 4.5 GeV) of 
the bottom quark, the mass (m^ = 175 GeV) of the top quark, the mass (m^o = 91.2 GeV) 
of the the mass {Ma = 200 - 400 GeV) of the CP-odd Higgs scalar and tan/3 = 35. 



A. Triple Higgs Bosons Self-Coupling Via e+e bbhh{HH, AA) 

To illustrate our results on the sensitivity to the hhh, Hhh, HHh, HHH, hAA, HAA 
triple modified Higgs bosons self-coupling, we show the n dependence of the total 
cross-section for e~^e~ bbhh{HH, AA) in Fig. |H We consider one representa- 

tive value of the Higgs boson mass. Ma = 400 GeV, and tan/? = 35 with a 
center-of-mass energy of i/i = 500, 1000, 1600 GeV and vary the triple couplings 
i^^hhh, i^^Hhh, i^^HHh, i^^HHH, i^^hAA, i^^HAA withiu the range K = —1 and +2. In all cases, 
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the cross-section is sensitive to the value of the triple couplings, as well as in the case of the 
process e^e~ —>■ bbhh, where for large values of Ma (the decoupling limit), the corresponding 
MSSM triple coupling approaches the SM triple coupling [50|. However, the cross-section 
and its sensibility to Xhhh, ^hHH decrease with increasing collider energy for e~^e~ bbhh, 
while for e+e^ tiHH{AA), the sensitivity to XhHH, Xhhh and Xhaa, Xhaa increases with 
rising collider energy. In particular, we can see in Fig. [5] for e~^e~ — > bbhh{HH,AA), the 
maximum cross-section (for M4 = 400 GeV and tan (3 = 35) is reached for y/s ~ 500 GeV 
and 1800 GeV respectively, which is consistent with Ref. 50|] in the case where h is identi- 
fied as the Higgs boson of the standard model. As an indicator of the order of magnitude, 
in Tables [Il lIIII we present the Higgs boson number of events (we have to multiplicate by 
the corresponding Branching Ratios to obtain the observable number of events) for sev- 
eral Higgs boson masses Ma, center-of-mass energy and k, values and for a luminosity of 
1000/6"^. If we consider the h bb decay for < 130 GeV, we have some opportunity 
to detect this process. In this region, the number of events is small but sufficient to detect 
e~^e~ bbhh — > bbbbbb. The BR{h bb) ~ 0.6 and the background for 6 b-jet are small. 

Finally, for completeness, in Fig. [6] we include a contour plot for the number of events of 
the studied processes as a function of ^/s and n with M4 = 400 GeV and tan jS = 35. These 
contour are obtained from Tables [ll lIIII Because in major parts of the MSSM parameter 
space the neutral Higgs bosons H and A are quite heavy, it is difficult to detect the processes 
e'^e~ bbH H (AA) when the relevant mechanism is e^e~ — > bbhh, as shown in the number 
of events given in Tables [Il lIIII 



B. Triple Higgs Bosons Self-Coupling Via e+e" tthh{HH,AA) 

As in the case of the processes e+e~ —>■ bbhh{HH, AA), Fig. [7] shows the variation of the 
cross-sections for e+e~ tthh{HH, AA) with nXhhh, i^Xuhh, i^Xnuh, hXhhh, i^Xhaa, i^Xraa, 
K = -1 to +2, 500 < v^s < 1600 GeV, Ma = 400 GeV and tan/3 = 35. The cross-section 
and its sensibility to the triple self-coupling on the energy range considered is dependent on 
the behavior of the cross-section with the center-of-mass energy determined by the position 
of the maximum. In Fig. [HI we can see that the maximum cross-section is reached for 
y/s ~ 1100 GeV in the case of e"'"e~ tihh, while for e"'"e~ tiHH{AA) it is reached for 
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TABLE I: Total production of Higgs boson pairs in the MSSM for tan /? = 35, £ = 1000 fb'^ and 
K = 0.5. 



Total Production of Higgs Boson Pairs 


e+e- ^ bbhh{HH, AA) k = 0.5 










MA{GeV) 


500 GeV 


1000 GeV 


1600 GeV 


200 


26 (-,-) 


17 (7,7) 


10 (7,7) 


250 


26 (-,-) 


17 (7,7) 


10 (7,7) 


300 


26 (-,-) 


17 (7,7) 


10 (7,7) 


350 


26 (-,-) 


17 (7,7) 


10 (7,7) 


400 


26 (-,-) 


17 (7,7) 


10 (7,7) 



TABLE II: Total production of Higgs boson pairs in the MSSM for tan (3 = 35, C = 1000 fb'^ and 
K = 1{MSSM). 



Total Production of Higgs Boson Pairs 


e+e- bbhh{HH, AA) n = 1{MSSM) 










MA{GeV) 


500 GeV 


1000 GeV 


1600 GeV 


200 


33 (-,-) 


19 (7,7) 


11 (7,7) 


250 


33 (-,-) 


19 (7,7) 


11 (7,7) 


300 


33 (-,-) 


19 (7,7) 


11 (7,7) 


350 


33 (-,-) 


19 (7,7) 


11 (7,7) 


400 


33 (-,-) 


19 (7,7) 


11 (7,7) 



^/s ~ 1800 GeV . For the process e+e" tihh (see Fig. [HI Ref. [50[) and for large values 
of Ma (the decoupling limit), the corresponding MSSM triple coupling approaches the SM 
triple coupling 50|. On the other hand, for ttHH{AA), the cross-section of the order 10~^ 
fb is small and is difficult to measure in the collider. In Tables HVllVIl to indicate the order 
of magnitude, we present the Higgs boson number of events (we have to multiplicate for the 
corresponding Branching Ratios to obtain the observable number of events) for several Higgs 
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TABLE III: Total production of Higgs boson pairs in the MSSM for tan/3 = 35, £ = 1000 /6^^ 
and K. = 1.5. 



Total Production of Higgs Boson Pairs 


e+e- ^ bbhh{HH, AA) k = 1.5 








V^ = 


MA{GeV) 


500 GeV 


1000 GeV 


1600 GeV 


200 


42 (-,-) 


22 (7,7) 


12 (7,7) 


250 


42 (-,-) 


22 (7,7) 


12 (7,7) 


300 


42 (-,-) 


22 (7,7) 


12 (7,7) 


350 


42 (-,-) 


22 (7,7) 


12 (7,7) 


400 


42 (-,-) 


22 (7,7) 


12 (7,7) 



boson masses Ma, center-of-mass energy and k values and for a luminosity of 1000/6"^. 

For e~^e~ — ^ tihh, the most favorable situation is for a center-of-mass energy of 1100 GeV 
and Mfi < 130 GeV, but in this case, the Branching Ratios for the four decay modes make 
this process very much suppressed. 

Finally, we include a contour plot for the number of events of the studied processes in the 
(-y/i, k) plane with Ma = 400 GeV and tan/5 = 35 in Fig. [9l These contours are obtained 
from Tables llVtlVIl Since in major parts of the MSSM parameter space the neutral Higgs 
bosons H and A are quite heavy, it is difficult to detect the processes e'^e~ — > ttHH{AA) 
when the relevant mechanism is e^e~ tthh for y/s ~ 1100 GeV as shown in Tables irvTlVIl 

C. Triple Higgs Boson Self-Coupling Via e+e^ bbhh{HH, AA),tthh{HH, AA) at 
CLIC energies 

In this subsection we analyze the triple Higgs self-coupling 
tiXhhh, i^^Hhh, f^^HHh, f^^HHH, i^'^hAA, i^^HAA via the proccsscs e~^e~ — > bbhh{HH, AA) and 
e^e~ — > tthh{HH,AA) for energies expected at the CLIC Figures fTOl and [TT] show the 
total cross-section for the double Higgs-strahlung in e~^e~ collisions, e~^e~ —>■ bbhh{HH, AA) 
and e~^e~ tthh as a function of k, for the cm. energy of y/s = 3 TeV, Ma = 400 GeV 
and tan /5 = 35. From these figures, we observe that the production cross-section of both 
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TABLE IV: Total production of Higgs boson pairs in the MSSM for tan/3 = 35, £ = 1000 fb~^ 
and K. = 0.5. 



Total Production of Higgs Boson Pairs 


e+e" tthh{HH, AA) k = 0.5 




V~s = 


V~s = 




MA{GeV) 


500 GeV 


1000 GeV 


1600 GeV 


200 




21 (-r) 


17 (-,-) 


250 




21 (-,-) 


17 (-,-) 


300 




21 (-r) 


17 (-,-) 


350 




21 (-,-) 


17 (-,-) 


400 




21 (-,-) 


17 (-,-) 



TABLE V: Total production of Higgs boson pairs in the MSSM for tan /? = 35, £ = 1000 fb''^ and 
K = 1{MSSM). 



Total Production of Higgs 


Boson Pairs 


e+e" ^ tthh{HH,AA) 


K = 1{MSSM) 






V~s = 






V~s = 


MA{GeV) 




500 GeV 


1000 GeV 




1600 GeV 


200 






23 (-,-) 




19 (-,-) 


250 






23 (-,-) 




19 (-,-) 


300 






23 (-,-) 




19 (-,-) 


350 






23 (-,-) 




19 (-,-) 


400 






23 (-,-) 




19 (-r) 



processes is small, as it is of the order of 10^'^ fh. 

Finally, in Tables IVIII and IVIIII we present the Higgs boson {hh, HH, AA) number of 
events for several k, values, luminosities of 1000 and 5000 fb~^ and center-of-mass energy 
= 3 TeV (we have to multiplicate by the corresponding Branching Ratios to obtain 
the observable number of events). It is evident from Figures [lO] and [TT] and Table IVIII 
that it would be difficult to obtain a clear signal of the processes e+e~ — > bbhh{HH,AA), 
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TABLE VI: Total production of Higgs boson pairs in the MSSM for tan/3 = 35, £ = 1000 fb~^ 
and K. = 1.5. 



Total Production of Higgs Boson Pairs 


e+e" tthh{HH, AA) k = 1.5 




V~s = 


V~s = 




MA{GeV) 


500 GeV 


1000 GeV 


1600 GeV 


200 




26 (-,-) 


21 (-,-) 


250 




26 (-,-) 


21 (-,-) 


300 




26 (-,-) 


21 (-,-) 


350 




26 (-,-) 


21 (-,-) 


400 




26 (-,-) 


21 (-,-) 



TABLE VII: Total production of Higgs boson pairs in the MSSM for tan (3 = 35, Ma = 400 GeV, 
^ = 3 TeV and C = 1000 fb-^. 



Total Production of Higgs Boson Pairs 



K 


e^e 
bbhh 


e~^e 
bbHH 


e+e 
bbAA 


e'^e — > 
tthh 


0.5 


5 


3 


3 


7 


1 


5 


3 


3 


8 


1.5 


5 


3 


3 
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Ma = 400 GeV, = 3 TeV 



e~^e tthh{HH, AA) at energies of a future linear collider such as CLIC, after having 
considered the background, except for ^/s = 3 TeV and very high luminosity {C = 5000 



fb ^), as shown in Table |VT 
double Higgs fusion process 



However, for the CLIC center-of-mass energy, the WW 



30 



31 



32 



33 



34( 1 . which increases with rising ^/s, can 



be exploited by large energies and luminosities and would be the preferred mechanism to 
measure the triple Higgs self-couplings. 
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TABLE VIII: Total production of Higgs boson pairs in the MSSM for tan/3 = 35, Ma = 400 GeV , 
^ = 3 TeV and C = 5000 fh-^. 



Total Production of Higgs Boson Pairs 





e^e — > 


e~^e — > 


e^e — > 


e'^e — > 


K 


bbhh 


bbHH 


bbAA 


tthh 


0.5 


25 


15 


14 


35 


1 


26 


16 


15 


41 


1.5 


28 


18 


17 


47 



= 400 GeF, 



^ = 3 TeV 



IV. CONCLUSIONS 

e"'"e~ linear colliders represent a possible opportunity for the triple Higgs boson self- 
coupling analysis. Therefore, we have analyzed the neutral Higgs bosons self-couplings with 
a complete set of the tree-level Feynman diagrams in the framework of the MSSM. The 
dependence of the triple Higgs bosons self-coupling Xhhhi ^Hhh, ^hAA, ^haa, ^hHH and 
^HHH on tan P and energy y/s were evaluated. 

The extended Higgs spectrum in supersymmetric theories gives rise to a plethora of 
triple couplings. The hhh coupling is generally quite different from the one expected in the 
standard model. It can be measured in Higgs double production at Future International 
Linear e~^e~ Colliders (ILC). Even though the e'^e~ cross sections are below the hadronic 
cross sections, the strongly reduced number of background events renders the search easier 
for the Higgs-pair signal, through bbbb final states for instance, in the e~^e~ environment 
than in jetty LHC final states. For sufficiently high luminosities the first phase of these 
colliders with an energy of 500 GeV will allow the experimental analysis of self-couplings 
for Higgs bosons in the intermediate mass range. Other couplings between heavy and light 
MSSM Higgs bosons can be measured as well, though only in restricted areas of the (v^, k) 
plane as illustrated in the sets of Figs. [6] and [9l 

In summary, we have analyzed the triple Higgs bosons self-coupling at future e~^e~ col- 
liders energies, with the reactions e~^e~ — bbhh{HH, AA) and e~^e~ —>■ tihh{HH, AA). The 
first is the most important after considering the h decay, and is statistically sufficient for an 
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accurate determination of k. In the case of e^e~ — > tthh{HH,AA), after we consider the 
t and h decays the final number of events is small; although they produce a non-vanishing 
number of events, statistically it is insufficient for the determination of k. In these condi- 
tions the number of events is small but our results have never been reported in the literature 
before and could be of relevance for the scientific community. 
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FIG. 1: Feynman diagrams at tree-level for e^e — > The diagrams for e^e — > tt/i/i are 

similar. 16 
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FIG. 2: Feynman diagrams at tree-level for e~^e bbHH. The diagrams for e^e ^ ttHH are 
similar. - ^ 
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FIG. 3: Feynman diagrams at tree- level for e^e — > hhAA. The diagrams for e^e — > ity4y4 are 
similar. 
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FIG. 4: Variation of the cross-section a{bbhh, bbHH, bbAA) with the modified triple coupling nXhhhj 
KXuhh at a colhder energy oi ^/s = 500, 1000, 1600 GeV with Ma = 400 GeV and tan/3 = 35. 
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FIG. 5: The dependence of the cross-section on center-of-mass energy y/s for Ma = 400 GeV 
and tan/3 = 35. The variation of the cross-section for modified triple couplings nXhhhj f^^Hhh is 
indicated by k = 0.5, 1.5. 
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function of ^/s and k. 
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25 




26 



